Detection of D2H + in the Dense Interstellar Medium 
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The 692 GHz para ground-state line of D2H + has been detected at the Cal- 
lO ■ tech Submillimeter Observatory towards the pre-stellar core 16293E. The derived 



ABSTRACT 



D 2 H + abundance is comparable to that of H 2 D + , as determined by observations 



CO 
O 

of the 372 GHz line of ortho-H2D + . This is an observational verification of recent 
theoretical predictions (Roberts, Herbst & Millar 2003), developed to explain the 
large deuteration ratios observed in cold, high-density regions of the interstellar 
medium associated with low mass pre-stellar cores and protostars. This detection 
confirms expectations that the multiply deuterated forms of H3 + were missing 
factors of earlier models. The inclusion of D 2 H + and D3 + in the models leads to 
predictions of higher values of the D/H ratio in the gas phase. 
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1. Introduction 

Recently, millimeter and submillimeter spectroscopy of the dense interstellar medium 
has shown that, in cold dense regions, deuterated molecular species are highly abundant, 
sometimes more than 10 _1 of the H version. Amazingly, doubly and triply deuterated species 
can be observed, e.g. D 2 CO (Ceccarelli et al. 1998), NHD 2 (Roueff et al. 2000), CHD 2 OH 
(Parise et al. 2002), D 2 S (Vastel et al. 2003), ND 3 (Lis et al. 2002; van der Tak et al. 2002), 
CD 3 OH (Parise et al. 2004). Several models have been developed to account for such high 
levels of deuteration (Tielens 1983; Roberts & Millar 2000a,b). Phillips & Vastel (2002) 
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have pointed out that the deuteration of H 3 + will be extended beyond H 2 D + , to D 2 H + and 
D 3 + , and that detection of D 2 H + might be possible. A calculation taking a high degree 
of deuteration into account has been carried out by Roberts, Herbst & Millar (2003) and 
Walmsley, Flower & Pineau des Forets (2004), confirming the expectation that, in dense 
depleted regions, the abundance of D 2 H + will be similar to that of H 2 D + , and that D3 + will 
be abundant. 

The key enabling work in the astronomical search for D 2 H + is the laboratory measure- 
ment of the para ground-state transition (l 10 -loi) by Hirao and Amano (2003). We report 
here the first astronomical detection of that transition. 

Chemical reactions go in the direction to minimize energy. The chemical fractionation 
process favors the production of the heavier more deuterated species, because of the mass 
dependence of the zero-point vibration energies of the isotopic variants. Gas phase species 
are expected to be depleted at the centers of cold, dark clouds, since they accrete on the dust 
grains (see, e.g., Charnley 1997). A series of observations has shown that the abundances of 
molecules like CO decrease in many pre-stellar cores (Bacmann et al. 2002). The removal 
of these reactive species affects the gas-phase chemistry and particularly the deuterium 
fractionation within the cloud. Indeed, the removal of species that would normally destroy 
H3 + (e.g. CO; Roberts & Millar 2000a) means that H3 + is more likely to react with HD and 
produce H 2 D + . For example, if [CO/H 2 ] ~ 5 x 10 -6 (Bacmann et al. 2002), this leaves HD 
at [HD/H 2 ] ~ 5 x 10 -5 as the most abundant molecule available for reaction with and 
H 2 D + , and favors the production of high deuterium content molecules: 



where AE a , AE& and AE C are the released energies of the exothermic reactions. Using 
the zero-point energies computed by Ramanlal, Polyanski & Tennyson (2003), and the energy 
of the first allowed rotational state of the Hjj" molecule permitted by the Pauli exclusion 
principle (~ 92 K), these values are: AE a = ~ 230 K, AE 6 = - 180 K and AE C = ~ 230 K. 

After a long frustrating search (Phillips et al. 1985; Pagani et al. 1992; van Dishoeck et 
al. 1992; Boreiko & Betz 1993), and with the advent of new submillimeter receivers, H 2 D + 
was detected toward two young stellar objects, NGC 1333 IRAS 4A (Stark et al. 1999) 
and IRAS 16293-2422 A (Stark et al. 2004), although with relatively low signal strength. 



H+ + HD< — ► H 2 D + + H 2 + AE a 
H 2 D + + HD< — > D 2 H + + H 2 + AE b 
D 2 H + + HD < — ► D+ + H 2 + AE C 



(1) 
(2) 
(3) 
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The H 2 D + search has now been extended to pre-stellar cores, and has been detected with 
relatively strong emission (Caselli et al. 2003; Caselli et al. 2004, in preparation; Vastel et 
al. 2004, in preparation) confirming that H 2 D + is dramatically enhanced in a gas depleted 
of most molecules. 

The ammonia and DCO + emission around the proto-binary system IRAS 16293-2422 
(Wootten & Loren 1987; Mizuno et al. 1990; Lis et al. 2002) does not only peak on IRAS16293- 
2422 itself but shows a second peak, about 90" to the southeast, in a condensation called 
16293E. Lis et al. (2002) found that CO in this region is depleted by a factor of 7. It is 
known to be a region where deuterium fractionation is strong and was chosen to be searched 
for D 2 H+. 

2. Observations and results 

The pure rotational transition (l 10 -loi) of D 2 H + has been measured in the laboratory 
by Hirao and Amano (2003). Spectroscopic observations of 16293E, presented here, were 
carried out in February 2004 using the facility receivers and spectrometers of the Caltech 
Submillimeter Observatory (CSO) on Mauna Kea, Hawaii. The position chosen was the 
DCO+ peak emission (a 200 o=16h32m29.4s, 5 20 oo=-24 28'52.6", Lis et al. 2002). We observed 
both the lio-ln transition of ortho-H 2 D+ (v = 372.42134(20) GHz, Bogey et al. 1984) and 
the lio-loi transition of para-D 2 H+ (u = 691.660440(19) GHz, Hirao and Amano 2003) (see 
Figures 1 and 2) which are the only lines currently available for these species. The data were 
taken under good weather conditions (225 GHz zenith opacity between 0.03 and 0.05). The 
CSO main beam efficiencies are ~ 60% for the 345 GHz receiver and ~ 40% for the 650 GHz 
receiver, determined from total power observations of Mars and Saturn. If the emission is 
extended compared to the beam size of CSO, as appears to be the case for DCO + (Lis et 
al. 2002), then the efficiency is about 70% at 372 GHz and 60% at 692 GHz. The FWHM 
beam size at 372 GHz is about 20", compared to ~ 11" at 692 GHz. Typical calibration 
uncertainties are ~ 24%. The pointing of the telescope was determined from observations 
of Jupiter, and was stable about ~ 2" (rms). We used both the 50 MHz and 500 MHz 
bandwidth acousto-optical facility spectrometers. The 500 MHz system was used to check 
system performance and calibration with the CO (6 — > 5) line, which is offset by ~ 82 km s _1 
from D 2 H + . There are no known lines of other interstellar molecules within 50 MHz, likely to 
emit from such a cold region. In Table 1 we report the frequency, the antenna temperature, 
the linewidth and the velocity relative to the local standard of rest, for the two lines. In 
Figure 3, we present 8.2 a and 4.4 a detections of H 2 D + and D 2 H + respectively, obtained in 
23 minutes and 103 minutes on-source integration time, respectively. 
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Since the para ground-state transition for D 2 H + is the only line available to existing 
telescope facilities, it is not possible to obtain confirmation of the identification. However, 
the situation is very different from molecule detection in hot core regions, such as OMC-1, 
where line confusion is rampant. At the 692 GHz transition frequency, any heavy molecule 
would need to be in a high J state, but the 10 K excitation temperature means that any 
such state cannot be occupied. Of course, CO (6 — > 5) is observable, but only at T a * = 2.3 
K. No U-lines have been observed in this source. The identification rests on the comparison 
of the deduced Vlsr for D 2 H + with that for other deuterated molecules, and also the line 
width and line strength. Note that the quoted uncertainty for the line frequency is about 
0.16 km s" 1 for H 2 D + and 0.008 km s" 1 for D 2 H + . The very slight difference in the observed 
Vi sr for these two lines (~ 0.2 km s _1 ) could be due, in part, to the uncertainty in the 
line frequencies, particularly H 2 D + . Also the accuracy of the astronomical measurement is 
limited by the resolution of the acousto-optic spectrometer at about 0.1 km s _1 . 



3. Discussion 

From the observed line strengths, given in column 3 of Table 1, we estimate the H 2 D + 
and D 2 H + column densities (see Table 2) for an excitation temperature T ex of 10 K, assuming 
a 25% calibration uncertainty (3 a). The column density is given by: 

N, °< = — aA - i J Tdv (4) 

where Q(T e:r ) is the partition function, Assuming LTE conditions, we can estimate the 
optical depth from the observed line intensity: 



T mb = [J u (T ex ) - J„{T bg )](l - e~ T ) (5) 

where J V (T) = [hv/k)/{e hu ^ kT — 1) is the radiation temperature of a blackbody at 
a temperature T, and T bg is the cosmic background temperature of 2.7 K. In the case of 
the H 2 D + transition, g u = 9, A ul = 1.04 10~ 4 s _1 , E ul = 17.9 K; in the case of the D 2 H+ 
transition, g„ = 9, A u i = 4.55 10~ 4 s _1 , E u i = 33.2 K. The derived column densities depend on 
the assumed value of the excitation temperature. Using NH 3 , Mizuno et al. (1990) estimate 
the gas temperature to be 12 K. Using D 2 CO line ratios, Loinard et al. (2001) obtained 
a rotational temperature of 8-10 K. Thus we quote, in Table 2, the values obtained for 
an excitation temperature of 10 K. Figure 4 presents the evolution of the ortho-H 2 D + and 
para-D 2 H + column densities as well as the para-D 2 H + /ortho-H 2 D + ratio, as a function of 
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temperature between 9 and 15 K. Figure 4 and Table 2 represent the case where the source 
emission is extended compared to the beam size. If the source emission is comparable to or 
smaller than the CSO beam size, the para-D 2 H + /ortho-H 2 D + ratio is then increased by a 
factor of ~ 1.5 at the average excitation temperature of 10 K. At thermal equilibrium, the 
ortho to para (respectively para to ortho) concentration ratio for H2D + (respectively D2H + ) 
is equal to 9 x exp(-86.4/T) (respectively 9/6 x exp(-50.2/T)), so that at 10 K, this ratio 
would be ~ 2 x 10~ 3 (respectively ~ 10~ 2 ). However, taking into account the limited rates 
of the spin allowed collisions with H 2 , it is found that at these low temperatures, the ortho to 
para H 2 D + concentration ratio is close to unity (Gerlich, Herbst & Roueff 2002). The para 
to ortho D 2 H + ratio is estimated by Walmsley, Flower & Pineau des Forets (2004) to be 
about one, for the same conditions. The para-D 2 H + /ortho-H 2 D + ratio presented in Figure 
4 should then approximately represent the actual D 2 H + /H 2 D + ratio. 

The 1.3 mm dust continuum strength (see Lis et al. 2002) is ~ 0.3 Jy in a 11" beam 
and ~ 1.3 Jy in a 20" beam corresponding to the angular resolution of the H 2 D + and D 2 H + 
data. Assuming a dust temperature of 12 K and a mass opacity coefficient of 0.005 cm 2 g _1 
(appropriate for pre-stellar cores), we derive an H 2 column density of ~ 5 1 23 cm" 2 . We 
then derive the H 2 D + and D 2 H + abundances to range between ~ 10~ 10 (at 10 K) and ~ 
10 -11 (at 15 K) compatible with abundances found by Roberts, Herbst & Millar (2003) for 
a cloud at 10 K and n(H 2 )=3 10 6 cm" 3 . 

The main result of this work is the detection of D 2 H + , with an abundance comparable 
to that of H 2 D + . This is a remarkable verification of recent theoretical predictions, aimed at 
explaining the large deuteration ratios observed in low mass pre-stellar cores and protostars. 
Two models, Roberts, Herbst & Millar (2003) and Walmsley, Flower & Pineau des Forets 
(2004), have recently considered the effect of including all possible deuterated isotopomers 
of H3 in the chemical networks, as suggested by Phillips & Vastel (2002). Roberts, Herbst 
& Millar (2003) studied the temporal evolution of a cold and dense cloud, and found that at 
late times (~ 10 4 yr), when CO is severely depleted in the gas phase (more than a factor of 
1000), the D 2 H + /H 2 D + ratio reaches unity. Walmsley, Flower & Pineau des Forets (2004) 
studied the evolution of gas, depleted in CO, as a function of gas density and grain size 
distribution. For densities larger than 10 6 cm -3 , they also found that D 3 + can be the most 
abundant ion and that the D 2 H + /H 2 D + ratio reaches unity (see their figure 2). Both models 
need an extreme CO depletion to account for such a ratio. As discussed in the Introduction, 
the measured CO depletion in 16293E is a factor 7 (Lis et al. 2002), rather than the extreme 
CO depletion needed by the models. However, the CO was measured in a large region (31") 
compared with that probed by the present observations: 20" for the H 2 D + and 11" for the 
D 2 H + observations. Also, regions containing CO will not contain much H 2 D + and D 2 H + 
and vice- versa, so the inevitable inhomogeneities in the region inhibit a clear result. 
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In summary, after some years of inconclusive results for theoretical models in under- 
standing the observed high deuteration ratios of doubly and triply deuterated molecules, 
the present observation seems to suggest that the basic process is now almost completely 
understood: the large deuteration is due to extreme CO depletion, and the factor that was 
previously missing in the models is the multiply deuterated forms of H^. This is quite an 
achievement, and one remaining step will be to verify that the last prediction, a significant 
abundance of D3 , is also fulfilled. At present, the lines detected here are the only ones 
available for H 2 D + and D 2 H+. Knowledge of the ortho to para ratios and abundances of 
H 2 D + and D 2 H + would be considerably improved if the ground state transitions of para- 
H 2 D+ (at 1370.15 GHz) and ortho-D 2 H+ (at 1476.60 GHz) were available. These lines could 
be detected from space telescopes such as Herschel with the Heterodyne Instrument for the 
Far-Infrared and also possibly from the stratospheric observatory SOFIA. However, D 3 + , 
like H 3 + has no permanent dipole moment. Therefore this molecule probably can only be 
detected in absorption in the near-infrared. 

The authors are grateful to John Pearson for providing an earlier spectroscopic compu- 
tation for D 2 H + and to Cecilia Ceccarelli for discussion. This research has been supported 
by NSF grant AST-9980846 to the CSO. 
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Fig. 1. — Diagram of the lowest energy levels of the D 2 H + molecule. 
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Fig. 2. — Diagram of the lowest energy levels of the H 2 D + molecule. 




Fig. 3. — Spectra of the ortho-H 2 D + lio-ln and para-D 2 H + lio-loi transitions towards 
16293E. 
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Fig. 4. — Column densities of the ortho-H2D + lio-ln (solid line) and para-D2H + lio-loi 
(dashed line) on the left side of the graph, and the para-D 2 H + /ortho-H 2 D + ratio (dotted 
line) on the right side, in the case where the source emission is extended compared to the 
beam size (see text). 
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Table 1. Results of Gaussian fits to the H 2 D + and D 2 H + spectra. 



Line 




Av 






(GHz) (K) 


(km s _1 ) 


(km s _1 ) 


H 2 D+ (lio-ln) 


372.42134(20) a 1.31 


0.36 ± 0.04 


3.55 ± 0.02 


D 2 H+ (lio-loi) 


691.660440(19) b 0.34 


0.29 ± 0.07 


3.76 ± 0.03 



a Measured frequency by Bogey et al. (1984). 
b Measured frequency by Hirao and Amano (2003). 
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Table 2. Ortho-H 2 D + and para-D 2 H + column densities at T ex = 10 K, in the case where 
the source emission is extended compared to the beam size (see text). 





T ex = 10 K 


t N(10 13 cm" 2 ) 


o-H 2 D+ 


0.74 1.73 ± 0.43 


p-D 2 H+ 


0.61 1.29 ± 0.32 


p-D 2 H+/o-H 2 D+ 


0.75 ± 0.37 



